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ABSTRACT 



Aims. In the approximation of linear dissipative magnetohydrodynamics (MHD) it can be shown that driven MHD waves in magnetic 
plasmas with high Reynolds number exhibit a near resonant behaviour if the frequency of the wave becomes equal to the local Alfven 
(or slow) frequency of a magnetic surface. This near resonant behaviour is confined to a thin region, known as the dissipative layer, 
which embraces the resonant magnetic surface. Although driven MHD waves have small dimensionless amplitude far away from the 
resonant surface, this near-resonant behaviour in the dissipative layer may cause a breakdown of linear theory. Our aim is to study the 
nonlinear effects in Alfven dissipative layer 

Methods. In the present paper, the method of simplified matched asymptotic expansions developed for nonlinear slow resonant waves 
is used to describe nonlinear effects inside the Alfven dissipative layer. 

Results. The nonlinear corrections to resonant waves in the Alfven dissipative layer are derived and it is proved that at the Alfven 
resonance (with isotropic/anisotropic dissipation) wave dynamics can be described by the linear theory with great accuracy. 

Key words. Magnetohydrodynamics (MHD) - Methods: analytical - Sun: atmosphere - Sun: oscillations 



1. Introduction 

I Magnetic fields are ubiquitously present in solar and space plas- 
mas. For regions where plasma-beta (the ratio of the kinetic and 
magnetic pressures) is less than one, magnetism controls the dy- 
namics, topology and thermal state of the plasma. The magnetic 
field in the solar atmosphere is not dispersed, but it tends to 

' accumulate in thinner or thicker entities often approximated as 
magnetic flux tubes. These magnetic flux tubes serve as an ideal 
medium for guided wave propagation. 

One particular aspect of the solar physics that has attracted 
much attention since the 1940s is the very high temperature of 
the solar corona compared with the much cooler lower regions 
of the solar atmosphere requesting the existence of some mech- 
anism(s) which keeps the solar corona hot against the radiative 
cooling. One of the possible theories proposed is the transfer 
of omnipresent waves' energy into thermal energy by resonant 
absorption or resonant cou pling of waves (see e.g. iPoedts et alj 
ri990l:ISakurai et alll99lHGoossens et al.lll995h . 

Waves which were initially obs erved sporadically r tiainly 
in radio wavelengths (see e.g., iKai & Takavanagil 119731; 
lAschwanden et al] Il992l) are now observed in abundance 



in all wa velengths, especially in 
e.g. , iDeForest & Gurman 1998; 



(extreme) ultraviolet (see 
Aschwanden et al. 1999; 
Nakariakov et al. 1999; Robbrecht et al. .2001;, King et aL,2003i: 
Erdelyi & Taroyan 2008; Mariska et an i2008h . Since the plasma 
is non-ideal, waves can lose their energy through transport 
processes, however, the time over which the waves dissipate 
their energy is far too long. In order to have an effective and 
localized energy conversion, the plasma must exhibit transversal 
inhomogeneities relative to the direction of the ambient mag- 
netic field. It was recognized a long time ago that solar and space 
plasmas are inhomogeneous, with physical properties varying 
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over length scales much smaller than the scales determined 
by the gravitational stratification. Homogenous plasmas have 
a spectrum of linear eigenmodes which can be divided into 
slow, fast and Alfven subspectra. The slow and fast subspectra 
have discrete eigenmodes whereas the Alfven subspectrum is 
infinitely degenerated. When an inhomogeneity is introduced 
the three subspectra are changed. The infinite degeneracy of 
the Alfven point spectrum is lifted and replaced by the Alfven 
continuum along with the possibility of discrete Alfven modes 
occurring, the accumulation point of the slow magnetoacoustic 
eigenvalues is spread out into the slow continuum and a number 
of discrete slow modes may occur, and the fast mag netoacoustic 
point spect rum accumulates at infinity (see e.g., iGoedbloedl 
I1975LI1984 . 

According to the accepted wave theories, effective energy 
transfer between an energy caiTying wave and the plasma oc- 
curs if the frequency of the wave matches one of the frequencies 
in the slow or Alfven continua, i.e. at the slow or Alfven reso- 
nances. The Alfven resonance has been more frequently associ- 
ated with heating of coronal structures given the low-yS regime of 
the solar corona. Nevertheless, slow resonance cannot be ruled 
out as an additional source of energy transfer. From a mathemat- 
ical point of view, a resonance is equivalent to regular singular 
points in the equations describing the dynamics of waves, but 
these singularities can be removed by, e.g., dissipation. Recently 
resonant absorption has acquired a new applicability when the 
observed damping of waves and oscillations in coronal loops 
has been attributed to resonant absorption. Hence, resonant ab- 
sorption has become a fundamental constituent block of one 
of the newest bra nches of solar physics, called coronal seis- 
mology (see e.g., Nakariakov et al. 1999; Ruderman & Roberts] 
2002; Goosse ns et al . 2002; Arregui et al. 2007; Ballai et al| 
20081; iGoossens et al.ll2008l: Pferadas et aI.H2008 ) when applied 
to corona and solar magneto-seismology when applied to the en- 



2 



Clack et al.: Nonlinear effects in Alfven dissipadve layers 



tire coupled solar atmos phere (see e.g. jErdelvi et al.l2007l : IVerflil 
l2007tlVerth etaPIIOOTh . 

Given the complexity of the mathematical approach, most 
theories describing resonant waves are limited to the linear 
regime. Perturbations, in these theories, are considered to be 
just small deviations from an equilibrium despite the highly non- 
linear character of MHD equations describing the dynamics of 
waves and the complicated interaction between waves and plas- 
mas. Initial numerical i nvestigations of resonant waves in a non- 
linear limit (see e.g., Ofman & Davila|[T99 5) unveiled that the 
account of nonlinearity introduces new physical effects which 
cannot be described in the linear framework. 

The first attempts to describe the nonlinea r resonant waves 
analytic ally appeared after the papers by iRuderman et al.l 
(Il997allbl) which were followed by furthe r analysis by, e.g .. 



Ballai et al.' ( 1998): Ba llai & Erdelvil (ll999h : lRuderm"Mil hoOOh : 
Clack & Ballai (2008); however, all these papers focused on the 
slow resonant waves only. These studies revealed that nonlinear- 
ity does affect the absorption of waves. In addition, the absorp- 
tion of wave momentum generates a mean shear flow which can 
influence the stability of resonant systems. 

The present paper is the first analytical study on the nonlin- 
ear resonant Alfven wave, where we obtain governing equations 
using techniques made familiar from previous studies on nonlin- 
ear slow resonant MHD waves. Before embarking on the actual 
derivation, let us carry out a qualitative discussion. First of all 
we should point out that in plasmas with high Reynolds num- 
bers (as in the solar corona) efficient dissipation only operates in 
a thin layer embracing the resonant surface. This layer is called 
the dissipation layer. This restriction on the effect of dissipation 
makes the problem more tractable from a mathematical point 
of view, as outside the dissipative layer the dynamics of waves 
is described by the ideal MHD. Dissipation is a key ingredient 
of the problem of resonance. As it was mentioned earlier dissi- 
pation removes singularities in mathematical solutions. From a 
physical point of view dissipation is important as it is the mech- 
anism which relaxes the accumulation of energy at the resonant 
surface and eventually contributes to the global process of heat- 
ing. 

It is important to stress that the choice of dissipation has 
to be related to the very physics which is described as differ- 
ent waves are sensitive to different dissipative mechanisms. Due 
to the dominant role of the magnetic field in the solar corona, 
transport processes are highly anisotropic. Possible dissipation 
mechanisms acting in coronal structures can be desc ribed within 
the framework of Braginskii's theory ('Braginskii' 1965) as it 
was shown in ap p Ucatio n s by. e.g. | Er delvi & Goossens ( 1995j); 
IOfman & Davila ( !l995h : iMocanu eTalJ t2008). Alfven waves 
are incompressible and transversal (in polarization), therefore, 
it is sensible to adopt shear viscosity and magnetic resistivity. 
Despite both transport processes being described by rather small 
coefficients, the net effect of dissipation can be increased consid- 
erably when the dissipative coefficients are multiplied by large 
transversal gradients. 

The paper is organized as follows. In the next section we 
introduce the fundamental equations and discuss the main as- 
sumptions. In Sect. Ill, we derive the governing equation for 
wave dynamics inside the Alfven dissipative layer Section IV 
is devoted to calculating the nonlinear corrections at Alfven res- 
onance. Finally, in Sect. V we summarize and draw our conclu- 
sions, pointing out a few applications and further studies to be 
carried out in the future. 



2. Fundamental equations and assumptions 

For describing mathematically the nonlinear resonant Alfven 
waves we use the visco-resistive MHD equations. In spite of 
the presence of dissipation we use the adiabatic equation as 
an approxim ation of the energy equation. Numerical studies by 
iPoedts et al.l ^1994) in linear MHD have shown that dissipation 
due to viscosity and finite electrical conductivity in the energy 
equation does not alter significantly the behaviour of resonant 
MHD waves in the driven problem. 

When the product of the ion (electron) gyrofrequency, 
and the ion (electron) collision time, t,((,), is much greater than 
one (as in the solar corona) the viscosity and finite electrical 
conductivity become anisotropic and viscosity is given by the 
Braginskii viscosity tensor (see Appendix A). The components 
of the viscosity tensor that remove the Alfven singularity are the 
shear components. The parallel and perpendicular components 
of anisotropic finite electrical conductivity only differ by a fac- 
tor of 2, therefore, we will consider only one of them without 
loss of generality. 

The dynamics of waves in our model is described by the 
visco-resistive MHD equations 



dp 

dt 



dt 



+ (v ■ V)v 



H- V ■ (pv) = 0, 

--VP+ — (B- V)B-H iv-S, 

P fJ-oP P 



5B 



— = V X (v X B) H- ^V^B, 

B2 

P = p + , V B = 0. 

2yU() 



(1) 

(2) 
(3) 
(4) 

(5) 



In Eqs. (dJl-® V and B are the velocity and magnetic induc- 
tion vectors, p the plasma pressure, p the density, A the coeffi- 
cient of magnetic diffusivity, y the adiabatic exponent, and po 
the magnetic permeability of free space. In addition, V ■ S is 
the Braginskii viscosity (see Appendix A for full details). Note 
that even though anisotropy has been considered the Hall term 
has been neglected. We neglect the Hall term from the induction 
equation, which can be of the order of diffusion term in the so- 
lar corona, because the largest Hall terms in the perpendicular 
direction relative to the ambient magnetic field identically can- 
cel. The components of the Hall term in the normal and parallel 
directions relative to the ambient magnetic field have no effect 
on the dynamics of Alfven waves in dissipative layers, hence 
these too are neglected. For full details on the Hall term and the 
reasoning behind neglecting it, we refer to Appendix B. 

We adopt Cartesian coordinates x, y, z and limit our analysis 
to a static background equilibrium (vo = 0). We assume that all 
equilibrium quantities depend on x only. The equilibrium mag- 
netic field. Bo, is unidirectional and lies in the yz-plane. The 
equilibrium quantities must satisfy the condition of total pres- 
sure balance. 



Po ■ 



- constant. 



(6) 



For simplicity we assume that the perturbations of all quan- 
tities are independent of y (d/dy = 0). We note that since 
the magnetic field is not aligned with the z-axis, Alfven waves 
still exist. In linear theory of driven waves all perturbed quan- 
tities oscillate with the same frequency, which means that 
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they can be Fourier-analyzed and taken to be proportional to 
exp(i[kz - cjt]). Solutions are sought in the form of propagating 
waves. All perturbations in these solutions depend on the combi- 
nation 9 = z-Vt, rather than z and t separately, with V - oj/k.ln 
order to match linear theory as closely as possible we apply the 
same procedure as above. In the context of resonant absorption 
the phase velocity, V, must match the projection of the Alfven 
velocity, va, onto the z-axis when x - xa where xa is the reso- 
nant position. To define the resonant position mathematically it 
is convenient to introduce the angle, a, between the z-axis and 
the direction of the equilibrium magnetic field, so that the com- 
ponents of the equilibrium magnetic field are 



Ov 



Bo sin a, B, 



Or 



Bo cos a. 



(7) 



The definition of the resonant position can now be written math 
ematically as 

V - Va (xA)cosa, 



where va is the Alfven speed defined as 

Bo 

VA 



In addition we introduce the sound and cusp speeds as 



7P() 
Po 



2 2 \l/2 
r V ' 



(8) 



(9) 



(10) 



In what follows we can take = without loss of generality. 
The perturbations of the physical quantities are defined by 



p^po+p, p ^ PQ+ p, B = Bo -H b, 

P 



Bo b b2 

P + + 



(11) 



where P is the perturbation of total pressure. 

The dominant dynamics of resonant Alfven waves, in lin- 
ear MHD, resides in the components of the perturbed mag- 
netic field and velocity that are perpendicular to the equilib- 
rium magnetic field and to the x-direction. This dominant be- 
haviour is created by an x"' singularity in the spatial solu- 
tion of these quantities at the Alfv en resonance (ISakurai et al.l 
119911: iGoossens & Rudermanlll995l) : these variables are known 
as large variables. The x- component of velocity, the compo- 
nents of magnetic field normal and parallel to the equilibrium 
magnetic field, plasma pressure and density are also singular, 
however, their singularity is proportional to ln|x|. In addition, 
the quantities P and the components of v and b that are parallel 
to the equilibrium magnetic field are regular; all these variables 
are called small variables. 

To make the mathematical analysis more concise and the 
physics more transparent we define the components of velocity 
and magnetic field that are in the yz-plane and are either parallel 
or perpendicular to the equilibrium magnetic field: 

















H 



w \l sma 
b, ] \ cos a 



V 

by 



w \ I COS a 
-b, ] \ sino- 



(12) 



where v, w, by and b- are the y- and z-components of the velocity 
and perturbation of magnetic field, respectively. 



Let us introduce the characteristic scale of inhomogeneity, 
/,„/,. The classical viscous Reynolds number, R^, and the mag- 
netic Reynolds number, are defined as 



_ PoVlinh Vlinh 

Kg - , Km - — = — , 

ri A 



(13) 



where po is a characteristic value of po, and ?; = ?7i is the shear 
viscosity coeflicient (see Appendix A). These two numbers de- 
termine the importance of viscosity and finite electrical conduc- 
tivity. We introduce the total Reynolds number as 



1 

R 



1 



1 

Rm 



(14) 



The aim of this paper is to study the nonlinear behaviour of 
driven Alfven resonant waves in the dissipative layer. We are not 
interested in MHD waves that have large amplitude everywhere 
and require a nonlinear description in the whole space. We fo- 
cus on waves that have small dimensionless amplitude e <K 1 far 
away from the ideal Alfven resonant point x = 0. 

In nonlinear theory, when studying resonant behaviour in the 
dissipative layer we must scale the dissipative coefficients (see 
e.g., Ruderman et al. 1997b; Ballai et al. 1998; Clack & Ballail 
I2OO8I) . The general scaling to be applied is 



^ = 7? '77, /I = 



(15) 



Linear theory predicts that the characteristic thickness of the dis- 
sipative layer, Idiss, is of the order of li„itR^^^^ and we assume 
that this is true in the nonlinear regime, too. Hence, we must in- 
troduce a stretching transversal coordinate, ^, in the dissipative 
layer defined as 



^^R 



1/3, 



We can rewrite Eqs. ([Tli-(|5]l in the scalar form as 



V 



dp d(pou) dw d(pu) d(pw) 



80 



dx 



•po 



de 



dx 



de 



(16) 



(17) 



du dP Bo cos a dbt _ I du dw \ 

PoV-- — I — — p u- — h w—— 

^ de dx po d0 ^\ dx do) 

du b_,db,, b;dbx _d^u 

~P^M ^ 'a^'^Tl' ^^^^ 

de p.Q dx pLo de dx^ 



d I Bocosff \ _/ dv_L dvA 

-^poVv,+Psma+— — ^,j=p^«— +w— j 

dv_i_ b_,db^ b.dbj_ _d^v^ 
de Po dx Po dO dx^ 



d I Bo cos a 

— po Vv|| - P cos ff + b\ 

m \ po 



dvi 



dvn 



bxdBo ^^^gv|| b,dbi\ b, db^ ^_d\ ^^^^ 
po dx de po dx po de dx^ ' 



Vbr + Bou cos a = wbr — ub, + A — - 

- ' de 



— (Vb^ + Bov± cos a) = — + — — — 

de dx de 



db,] 
dx J 



(21) 



_^^^_^^^_IV2^^, (22) 
dx ^ de 
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^ nzu , 73 ^ d{BQu) dw 
— [ybw + B(^v\\ cos a) -So — 

d{wb\i) dv\\ 5y|| - 2, 

= — 7. 1 T. b^— b-,— AN b\\, 

dx de dx ' de " 



(23) 



y 



'dp^ 

,de 



.dp 

' de 
1 

Po 



(dpo 


2 dpo 


\ dx 


'^^ dx 


^dp 


dp\ 
f'Tej- 


'de 



—dp _dp 

^^'de'^'de 



+u 



dpo 
dx 



_ dpo 
dx 



-dp 

dx 



5^ 

dx 



l-[bl + bl + bl + 2B^b^), 



db^ ^ dbz 

dx de 



(24) 

(25) 
(26) 



In the above equations V - {d/dx,0,d/de) and w - V||Cosq' - 
Vj_ sin a. 

Equations (fT7]i-(l26]l will be used in the following sections 
to derive the governing equation for the resonant Alfven waves 
inside the dissipative layer and to find the nonlinear corrections. 



3. The governing equation in the dissipative layer 

In order to derive the governing equation for wave motions in 
the Alfven dissipative l ayer we employ the method of rnatche d 
asymptotic expansions (lNavfehll98lHBender & Orszaghll991h . 
This method requires to find the so-called outer and inner expan- 
sions and then match them in the overlap regions. This nomen- 
clature is ideal for our situation. The outer expansion corre- 
sponds to the solution outside the dissipative layer and the in- 
ner expansion corresponds to the solution inside the dissipative 
layer. A simplified version of the method of matched asymptotic 
expansions, developed by Ballai et al] (1199 8), is adopted here. 

The typical largest quadratic nonlinear term in the system of 
MHD equations is of the form gdg/dz while the typical dissipa- 
tive term is of the form 7id^g/dz^, where g is any 'large' vaiiable. 
Linear theory predicts that 'large' variables have an ideal singu- 
larity jc"' in the vicinity of x - Q. This implies that the 'large' 
variables have dimensionless amplitudes in the dissipative layer 
of the order of eR^^^. It is now straightforward to estimate the 
ratio of a typical quadratic nonlinear and dissipative term. 



^ ^ gdg/dz 
* rjd^g/dz^ 



(27) 



where the quantity (f)^ can be considered as the quadratic nonlin- 
earity parameter. If the condition eR^^^ «: 1 is satisfied, linear 
theory is applicable. On the other hand, if eR^^^ > 1 then nonlin- 
earity has to be taken into account when studyin g resonant waves 
in dissi pative layers. Using the same scalings, iRuderman et alj 
( Il997bl) showed that nonlinearity has to be considered when- 
ever slow resonant waves are studied in the solar photosphere. 
For a typical dimensionless amplitude of e ~ 10"^ linear the- 
ory can be applied if the total Reynolds number is less than 10^. 
This value is much less than the resistive and shear viscosity 
Reynolds number (10'" - 10'^). This conclusion implies that 
in the solar atmosphere resonant absorption should be a non- 
linear phenomenon. In order to describe the role of dissipation 
and nonlinearity equally we assume that (pq ~ \. 



Far away from the dissipative layer the amplitudes of pertur- 
bations are small, so we use linear ideal MHD equations in order 
to describe the wave motion. The full set of nonlinear dissipative 
MHD equations are used for describing wave motion inside the 
dissipative layer where the amplitudes can be large. We, there- 
fore, look for solutions in the form of asymptotic expansions. 
The equilibrium quantities change only slightly across the dis- 
sipative layer so it is possible to approximate them by the first 
non-vanishing term in their Taylor series expansion with respect 
to X. Similar to linear theory, we assume that the expansions of 
equilibrium quantities are valid in a region embracing the ideal 
resonant position, which is assumed to be much wider than the 
dissipative layer This implies that there are two overlap regions, 
one to the left and one to the right of the dissipative layer, where 
both the outer (the solution to the linear ideal MHD equations) 
and inner (the solution to the nonlinear dissipative MHD equa- 
tions) solutions are valid. Hence, both solutions must coincide in 
the overlap regions which provides the matching conditions. 

Before deriving the nonlinear governing equation we ought 
to make a note. In linear theory, perturbations of physical quan- 
tities are harmonic functions of e and their mean values over a 
period are zero. In nonlinear theory, however, the perturbations 
of variables can have non-zero mean values as a result of non- 
linear interaction of different harmonics. Due to the absorption 
of wave momentum, a mean shear flow is generated outside the 
dissipative layer (Ofman & Davila 1995). This result is true for 
our analysis also, however, due to the length of this study we 
prefer to deal with this problem in a forthcoming paper 

We suppose that nonlinearity and dissipation are of the same 
order so we have eR^^^ - i.e. R ~ e^^^^. We can, therefore, 
substitute e^^^^ for R in Eq. ( fTSl l to rescale viscosity and finite 
electrical resistivity as 



,3/2 



A. 



(28) 



We do not rewrite the MHD equations as they are easily obtained 
from Eqs. ([T7j-(|26]l by substitution of Eq. jjsl ). 

The first step in our description is the derivation of govern- 
ing equations outside the dissipative layer where the dynamics 
is described by ideal {rf - A - Q) and linear MHD. The linear 
form of Eqs. (fT7]i-(|26]l can be obtained by assuming a regular 
expansion of variables of the form 

f^efUe'l^f^K.., (29) 

and collect only terms proportional to the small parameter e. 
This leads to a system of linear equations for the variables with 
superscript '1'. All variables can be eliminated in favour of m*'' 
and f leading to the system 



V- 



de 



dx 



dP'-'^ du^'^ 

V—:— =poA- 



where 



dx 



poAC 



de 



y4 - y2 (v^ +cfj + v]cl cos2 a ' 
A-V^-v\ cos^ a, 
C = [v\+cl){y^ -c\ cosV). 



(30) 



(31) 



(32) 



The quantities A and C vanish at the Alfven and slow resonant 
positions, respectively. As a result these two positions are regular 
singular points for the system ( l30b . The remaining variables can 
be expressed in terms of m^^' and P^'^ as. 



„(i) 



y sina^ 
PoA 



>(i) 



„(i) 



Vc\ cos a ,, , 
PoC 



(33) 



Clack et al.: Nonlinear effects in Alfven dissipative layers 



5 



80 pi)C 

dpW y2 „(i) 



(34) 



50 (35) 



(36) 



(37) 



561 C de V dx 

Since Eq. (|30] l has regular singular points, the solutions can be 
obtained in te rms of Frobenius series with respect to x ( for de- 
tails see, e.g.. iRuderman et alJ[T997bt iBallai et al.lll998l) of the 
form 



of the order of e In e in the dissipative layer Strictly speaking, 
the inner expansions of all variables have to contain terms pro- 
portional to eln e and e^'^ In e (see e.g., Ruderman et al. 1997b!). 
In the simplified version of matched asymptotic expansions we 
utilize the fact that | In e| <K e^" for any positive k a nd e — > -i-O, 
and c onsider In e as a quantity of the order of unity dBallai et al.l 
'1998). This enables us to write the inner expansions for u, b^, b\\, 
p and p in the form of Eq. ( |29] ). 

We now substitute the expansion ( |29] l for P, u, b^, b\\, vy, p 
and p and the expansion given by Eq. ( l42l i for and b_i_ into 
the set of equations obtained from Eqs. (fT7])-(|26b after substitu- 
tion of X = e'^^^. The first order approximation (terms propor- 
tional to e), yields a linear homogeneous system of equations for 
the terms with superscript ' 1' . The important result that follows 
from this set of equations is that 



p(i) ^ + P'^>(6)xln \x\ + P\"i0) + . . 



d(1)/ 



d(1)/ 



3(1)/ 



= m'i"(6') In + u'^'(6) + u'^"(6)xln\x\ + . 



,.(1) 



(1)/ 



,.(1)/ 



(38) 
(39) 



The coefficient functions depending on in the above expansions 
are, generally, different for x < and x > 0. The particular form 
of these series solutions indicates that the perturbation of the 
total pressure is regular at the ideal resonant position. From Eqs. 
(I33ll-(l37ll. we see that the quantity v||'' is also regular, while all 

other quantities are singular The quantities foj'', 

and p*^'* behave as In |x|, while v^^^ and /7*|' behave as x"', so they 
are the most singular 

As the characteristic scale of dissipation is of the order of 
linhR~^^^ and we have assumed that R ~ e^^^ we obtain that the 
thickness of the dissipative layer is li„i,R^^^^ = ^(e'^^Z,-,,/,), im- 
plying the introduction of a new stretched variable to replace the 
transversal coordinate in the dissipative layer, which is defined 
as ^ = e '^-x. Again, for brevity, Eqs. (fT7Ti-(l26]l are not rewritten 
as they can be obtained by the substitution of 



d_ 
dx 



.-1/2 



d_ 



(40) 



for all derivatives. The equilibrium quantities still depend on x, 
not ^ (their expression is valid in a wider region than the charac- 
teristic thickness of the dissipative layer). All equilibrium quan- 
tities are expanded around the ideal resonant position, x = 0, 
as 



dx 



(41) 



where fo is any equilibrium quantity and the subscript 'A' in- 
dicates that the equilibrium quantity is evaluated at the Alfven 
resonant point. 

We seek the solution to the set of equations obtained from 
Eqs. (fT7Ti-(|26]l by the substitution of x = e'''^^ into variables in 
the form of power series of e. These equations contain powers 
of e'^^, so we use this quantity as an expansion parameter To 
derive the form of the inner expansions of different quantities 
we have to analyze the outer solutions. First, since vy and P are 
regular at x = we can write their inner expansions in the form 
of their outer expansions Eq. ( |29] l. The amplitudes of large vari- 
ables in the dissipative layer are of the order of e'^^, so the inner 
expansion of the variables Vj^ and b± is 



p(i) ^ 



(43) 



that is P^^^ does not change across the dissipati ve layer. This 
result parallels the result fou nd in linear theory jSakurai et al.l 
119911 : iGoossens et alj [T995h and nonUnea r theories of slovy 
resonance (see e.g.. 'Ruderman e t al1ll997bl: iBallai et alJ[T998l: 
Clack & Ballai 2008). Subsequently, all remaining variables can 
be expressed in terms of u^^^ 



v|l' and as 



„(i) 



Cr COS a 
— 1 

vl poV 



>(i) 



V 



(44) 



(45) 



_Bo{vl-cl)dpm_^^ idBA 
d0 V \ dxj' 



30 



80 



v\ d0 



V fio \ dx 



5p<i> _ 1 m(i> Idpo 

~d0~ ~ vl d0 ^~\~dx 



(47) 



(48) 



All equilibrium quantities are calculated at x = 0. In addition, the 
relation that connects the normal and perpendicular components 
of velocity is 

^ - (49, 



■ sma- 



80 



0. 



In the second order approximation we only use the expres- 
sions obtained from Eqs. ( fT9] l and ( |22] |. Employing Eqs. ( |43] |- 
(|49] |. we replace the variables in the second order approximation 
which have superscript ' 1 ' . The equations obtained in the second 
order are 



8P'-^^ . Bo cos a 8b? 8v? 
■ sin Q- H + Vpo - 



80 



80 

„(i) 



80 



pav\ \ dx 80 \ dx 80 8f 



2.,(1) 



(50) 



^1/2^(1) + eg(2) + , 



(42) 



The quantities u, bx, b\\, p and p behave as In |x| in the vicinity 
of X = 0, which suggests that they have expansions with terms 



d_B,\ 



8v' 



d) 



|2„(1) 



8bf 8vf 

^^r^+^ocosa— ^-Hcosa \f—^-A^ 

80 80 \ dx 80 vUosa 8f 



BoV 8Wl 



(51) 
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Once the variables with superscript '2' have been eliminated 
from the above two equations, the governing equation for res- 
onant Alfven waves inside the dissipative layer is derived as 



V sin a dP'^^ 



Po 



de 



where 



^dvl-^ 



dx ^ 



cos^ a. 



(52) 



(53) 



It is clear that Eq. (1521 ) does not contain nonlinear terms de- 
spite considering the full MHD system of equations. This re- 
sult is in stark contrast with the results obtained for nonlinear 
slow resonance where the governing equation was found to be 
always nonlinear (see e.g., iRuderman etanil997bt[Ballai etal.1 
Il998i IClack & Baliaill2008h . The governing equation ^ sug- 
gests that resonant Alfven waves can be described by the linear 
theory unless their amplitudes inside the dissipative layer is of 
the order of unity. 

As the quadratic nonlinear terms cancel each other out, it 
is natural to take into account cubic nonlinearity (the system of 
MHD equations contain cubic nonhnear terms), where the non- 
linearity parameter is defined as 



g^dgldz 2r, 

6 K. 



rjd^gldx 



(54) 



Despite the higher order nonlinearity the governing equa- 
tion is similar to the equation derived for quadratic nonlinearity 
( [52] l. These results require finding an explanation to the linear 
behaviour of waves inside the dissipative layer. The following 
section will be devoted to the study of nonhnear corrections in 
the Alfven dissipative layer. 



4. Nonlinear corrections in the Alfven dissipative 
layer 

Since we have assumed that waves have small dimensionless 
amplitude outside the dissipative layer, we will concentrate only 
on the solutions inside the dissipative layer. 

In our analysis we use the assumptions and equations pre- 
sented in Section II, however, we will not impose any relation 
between e and R. Equations JTSl l and (fT&t will be used to de- 
fine the scaled dissipative coefficients and stretching transver- 
sal coordinate in the dissipative layer For simplicity we denote 
6 - R^^l^ . This means that our scaled dissipative coefficients and 
stretched transversal coordinate become 



(55) 



The first step to accomplish our task is to rewrite Eqs. (fT7]i-(l26]l 
by substituting 



d , d d d d d 

^^S'^, ^ = and — = -y— . 

dx d£, dz de dt de 



(56) 



All equilibrium quantities (which are still dependent on x, not 
^) will be approximated by the first non- vanishing term of their 
Taylor expansion (see, Eq. (14111). 

The substitution of Eqs. (I4ni55ll56] l will transform Eqs. (fTTT l- 
( |26] l into 

du dpo d(pu) d 



dP b,db, 5 db, 

-rr (Bo cos a + /?:) — 

d^ juo d^ juo de 



^ du du ,77 d^u 
5{V -w) — -M— + 5 — 

' de d^ podf 



(58) 



dP . b,db^ 6 ,db^ 

6—sma+ — + — (Bo cos a -H Z?^) — - 

de fio d^ fj.() de 

^^S(V-w)— + u—-6^ — — , (59) 
de d^ po d^^ 



dP b, dbu 6 , , db« 

6— cos a (Bo cos a + b~) — — 

de jUo d^ fio " de 



6 dBo^ 
Ho dx 



6(V -w)—^-u—^+5^ 

de d^ po df 



(60) 



N , .du 

6(V -w) -— + 6 (Bo cos a + b:,) — 

de de 

52 



+ 6'A{^+S'^\b,^0, (61) 



5^2 



^dbj_ db^ (du di'ii \ dv_L 



dv^ 



+ 6(Bo + 77||) -±cosa + 6'A\^+6'^\b^^0, (62) 



52 



df 



^ 5^711 (dbn dBo\ dvii dv\\ 
de \d^ dx I d^ de 



■ (Bo + fell) 



du dvj_ . 1 
o— — sm a H- A 

d^ de 



d d 

6(y-w) — -u— 
de d^ 



5^ 



56*2; 



0, 

(63) 
(64) 

(65) 



The only condition we need to impose when deriving the 
nonlinear corrections to resonant Alfven waves in the dissipa- 
tive layer is imported from the linear theory which predicts that 
in the dissipative layer 'large' variables have dimensionless am- 
plitude of the order of e/?'^^ (see Section III). We assume that the 
dimensionless amplitudes of the linear approximation of 'large' 
variables (v±_ and bi_) in the dissipative layer are small, so that 



e«/;-i/3(^5). 



(66) 



This condition ensures that the oscillation amplitude remains 
small inside the dissipative layer. From a naive point of view 
the linear theory is applicable as soon as the oscillation ampli- 
tude is small. The example of slow resonant waves clearly shows 
that this is not the case. The nonlinear effects become important 
in the slow dissipative layer as soon as e ~ R^^^^, i.e. as soon 
as the oscillation amplitude in the dissipative layer, which is of 
the order of e/?''^, is of the order of '^^ <K 1. For example, 
in the corona perturbations with dimensionless amplitudes less 
than 10""* can be considered by this theory. From Eq. ( |27l ) we 
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would expect to see quadratic nonlinearity appear for waves with 
dimensionless amplitudes larger than 10 and from Eq. ( |54] i we 
would expect to see cubic nonlinearity appear for waves with di- 
mensionless amplitudes larger than 10"^. If we take e ^ R^^^^ 
we find that inside the dissipative layer we have dimensionless 
amplitudes of the order of unity. This causes a breakdown in 
our theory, and therefore another approach would have to be 
adopted. At this time we do not know of an analytical study 
which can carry out this task without considering the full non- 
linear MHD equations throughout the domain. 

We now assume that all perturbations can be written as a 
regular asymptotic expansion of the form 



(67) 



where /o(x) represents the equilibrium value. Substitution of ex- 
pansion (|67| i into the system (|57]i-(|65ll leads to a system of equa- 
tions which contains the small parameter 6. This observation in- 
spires us to look for the solution in the form of expansions with 
respect to S. In order to cast large and small variables in this de- 
scription we are going to use the following expansion for small 
variables (u, bx, v\\, b\\, p, p and P) 



(68) 



while large variables (vj^ and b^^ will be expanded according to 



■(2) 



/ii - S h, + h, + 



(69) 



The bar notation is used here to distinguish between these ex- 
pansions and the expansions used in the previous section. From 
this point on we drop the bar notation. 

Substituting Eqs. (|68]|-(|69]| into the system (l57])-(l65b. taking 
terms proportional to e and then only retaining terms with the 
lowest power of 6, results in the set of linear equations 



a (1) 

Po^TT sina -po 



du 



(1) 



de 



= 0, 



dF 



.(1) 



0, 



PoV- 



d) 



Bo cos a db^^l 



de 



0, 



dP' 



(I) 



86 



cosa-poV-^ 



80 



£l_ 



dx 



80 
Bo cos a 



(70) 



(71) 



(72) 



(1) 



80 



V- 



8b' 



(1) 



80 



+ Bo cosa- 



8u 



(1) 



y- 



8y 



(1) 



86 



+ Bq cos a- 



86 

8v'!^ 



80 



= 0, 



8u 



(1) 



PoV 



80 



■PoclV 



8^ 
86 



3(1) 



- Bo sin a- 



8v' 



(1) 
±1 



86 



= 0, 



dpo 
dx 



dpo 
dx 



0. 



= 0, (73) 
(74) 
(75) 
(76) 

= 0, (77) 
(78) 



In these equations all equilibrium quantities are calculated at x = 
0. 



Using these equations we can express all dependent variables 



in terms of m*'', v^'J and 



(1) _ c^cosa (1) 



vl PoV 



(1) _ BqV 



Bocoso' n) 
V 1 



(79) 
(80) 



8b^ _Bo{.l-cl)dP\^^ ^u^l^idjo] ^3^^ 



80 



pov\ 



d6 V \ dx 



8p\'^ cj dP['^ Bo I dBo 



80 vl d0 V Ho \dx j' 



df^_]_dP^ uf^ (dpo 



It follows from Eq. dTTT l that 

Pf' = Pf\6). 
Finally, we obtain the relation between m*'* and v^'J, 



^ (1) a (1) 

8u\ 8v\[ 
~df ~ If 



sma. 



(82) 
(83) 

(84) 
(85) 



Note that Eqs. (TTgll-dSS]) ai-e formally identical to Eqs. ( |43] )- 
(|49] l for the linear approximation in Section III. This is not sur- 
prising as both methods are designed to replicate linear theory in 
the first order approximation. 

4. 1 . The second and third order nonlinear corrections 

Once the first order terms are known we can proceed to derive 
the second and third order approximations with respect to e (i.e. 
terms from the expansion of Eqs. (|57])-(|65]| that are proportional 
to and e^, respectively). First, we write out the second or- 
der approximations and substitute for all first order terms (i.e. 
terms of the form /['*) using Eqs. (|79]l-(l85]l. Secondly, we find 
(by solving the inhomogeneous system) the expansions of sec- 
ond order terms (terms with subscript '2'). Thirdly, we derive 
the second order relations between all variables, similar to the 
ones obtained in the first order approximation. 

The equations representing the second order approximation 
with respect to e (with variables in the first order substituted) are 



du2 ^ 
Po 1: — I- 



dpo 
dx 



8u2 



dpi I8v\\2 
V-— + Po —— cos a 
80 ^ \ 80 



8vj_2 . \ (dpo 

Sm Q- M- I M2 

86 \dx' 



v^dP^ 
'v\ d6 



sma ■ 



8P2 
8^ 



-6 



8u2 Bo cos a dbx2 

PoV— + — 

86 Ho 80 



(86) 
(87) 



8P2 8v_i_2 
— — smff + PoV— — 
561 ^ 86 



Bo cos a db^2 



HO 



86 



+ 6U 



dpo\ dvj^ 
dx] 80 



(dBo' 


8b_i_2 


\ dx 1 


86 



+ 77- 



5Vl2 

8f 



= ^(r'), (88) 
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— — cos a - pnV — — 

do ^ 36 



bx2 I dBp 
fi \ dx 

dpo \ dv\\2 



Bo cos a db\ 



l|2 



cos a I dB() \ db\\2 



dx I do 



cos a sin Q- ^\)dP^ 
V ^'^^~d9 



+ mi), (89) 



y- 



. dbx2 
86 



+ Bi\ cos a 



+ 6 



du2 

'de 

d^b,2 



+ ^cos a 





(9m 2 


m 


'do 



= ^(1), (90) 



5^7x2 dvjiji 
V— — I- Bq cos a- 



de 



de 



+ 6 



df 



+ ^cos a 





dv^2 


\ dx 1 


de 



dBo 
dx 



U2 - Bo sin a 



du2 
Bo^.5 



dv±2 

de 



V- 



dh 



I dBo' 


du2 


\ dx 1 


'di 



l|2 

de 

(I) 



Bo sin a dP^ 

povl 



de 



(92) 



P2- P2 



Mo 



-b\\2 + 6 



dBo 
dx 



dpo 
dx 

b\\2 



dpo 
dx 



= mSl, (93) 



PO (1) 



±1 



-'). 

(94) 

It is clear that nonlinear terms appear from this order of approx- 
imation and they are expressed in terms of variables obtained in 
the first order. 

The analysis of the system of Eqs. (I86]l-(l94li reveals that the 
expansions with respect to S has to be written in the form 



82 



for U2, bx2, v±2, b_i_2 and f 2 and 



h2^6-^h\'' + 6-'h?+h^'' + . 



(95) 



(96) 



fory||2, bp,, p2 andp2- 

Here we need to make a note. It follows from Eqs. ( |95] l 
and ( |96] l that the ratio of p2 to pi is of the order of e6^^, and 
the same is true for vp, bp and p2- It seems to be inconsistent 
with the regular perturbation method where it is assumed that 
the next order approximation is always smaller than the pre- 
vious one. However, this problem is only apparent. To show 
this we need to clarify the exact mathematical meaning of the 
statement "in the asymptotic expansion each subsequent term is 
much smaller than the previous one". To do this we introduce 
the nine-dimensional vector U - (m, vy, Vj., bx, b\\, b±, P, p,p) and 
consider it as an element of a Banach space. The norm in this 
space can be introduced in different ways. One possibility is 



Jo J-a 



|U| d^, 



(97) 



where L is the period. The asymptotic expansion in the dissipa- 
tive layer, Eq. ( l68T l, can be rewritten as U = Uo-HeUi -i-e^U2 + 



Then the mathematical formulation of the statement "each sub- 
sequent term is much smaller than the previous one" is | |U„+ 1 1 1 « 
||U„||, « = 1,2, It is straightforward to verify that, in accor- 
dance with Eqs. dMll, dSl, dSl and IIU2II « ||Ui||. 

Once the expansions ( |95] ) and ( |96] l are substituted into Eqs. 
(I86]l-(l94ll, we can express the variables in this order of approxi- 
mation as 

Bo - 



^\\2 



2V 



"112 



(1)^ 

2vf- 



^±2 



±2 



For the total pressure we obtain that 



dP^ 
d^ 



3(1) 



Bo cos a 



p'^\e). 



,(1) 



(98) 

(99) 
(100) 

(101) 

(1) 



') (91) addition, we obtain that the equation which determines is 



du 



(1) 



d^ 



V 



(102) 



Since the large variables in this order of approximation are 
V||2 and bp, we can deduce that the linear order of approxima- 
tion of resonant Alfven waves in the dissipative layer excite 
magnetoacoustic modes in the second order of approximation. 
The excitation comes from the nonlinear term found in the sec- 
ond order approximation of the pressure equation, this drives the 
parallel components of the velocity and magnetic field perturba- 
tions. Since we are focussed on the Alfven resonance only, these 
waves are not resonant. These waves act to cancel the very small 
pressure and density perturbations created by the first order ap- 
proximation. 

We now calculate the third order approximation with respect 
to e. On analyzing the third order system of equations we de- 
duce that the large variables in this order of approximation are 
Alfvenic {v% b^[!^), so we only need the perpendicular compo- 
nents of momentum and induction equations given by Eqs. ( |59l ) 
and ( |62] i, respectively. 

First, since some of the first order approximation terms con- 
tribute to the third order approximation in integral form we must 
introduce a new notation 



U 



(1) 



J u--^de. 



(103) 



The third order approximation of the perpendicular component 
of momentum is 



5P3 5vx3 Bo cos a db^^ 
— — sma -I- poV— 1 — h 

de ^" de po de 



\ \ dx j de po \ dx j 



db^3\ dh'^3 



= 6 



d^ 



u^'U?^ I dpo 



V 



dx 



de 
'), (104) 



while the perpendicular component of magnetic induction equa- 
tion is 



5/7 J.3 dv^ji 
V^r— + Bo cos Q-^r— + 6 



de 



de 



d^b^, 




de 




idBo\ 


5vj.3 


\ dx j 


de 



\ (105) 
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Note that in obtaining the third order approximations we have 
employed all the relations we have for variables in the first and 
second order of approximation. 

Equations jlOAj and dlOSt clearly show that the nonlinear 
terms on the right-hand sides do not cancel. This implies that the 
expansion of Vj_3 and i>j_3 should be of the form 

h3 = (J"^/!^" + d-^hf + + .... (106) 

We should state, for completeness, that if we derive the third 
order approximation for all the Eqs. (|57]i-(|65ll we obtain the ex- 
pansions for M3, bx-i, V||3, b\\},, p3, p3 and P3 to be 

8,^6-%'^^6-^8fUgf + .... (107) 

The expansions calculated for all the variables can now be 
collected together and we can write the expansions for 'large' 
and 'small' variables in the dissipative layer when studying res- 
onant Alfven waves. Large variables {vj_ and b±) have the ex- 
pansion 

= e(^)/.<'>.g)%<'>...., (108) 

and the expansion of small variables (m, b^, vy, b\\, p, p and P) is 
defined as 

g = .^<;)+(^)^^" + .(f)%<''H-.... (109) 

Since Eq. ( |66] ) is the only condition enforced in the dissipa- 
tive layer, we can state that 

'>(i)>(ir>(i)'>-- 

Therefore, since both Eqs. ( 1108b and ( |109l l contain succes- 
sive higher powers of the parameter e/6 we can deduce that, con- 
sidering Eq. ( II 10b . higher orders of approximation of large and 
small variables become increasingly insignificant in comparison 
to the linear order of approximation, so resonant Alfven waves in 
the dissipative layer can be described accurately by linear theory 
if condition ( |66] | is satisfied. 

5. Conclusions 

In the present paper we have investigated the nonlinear be- 
haviour of resonant Alfven waves in the dissipative layer in one- 
dimensional planar geometry in plasmas with anisotropic dissi- 
pative coefficients, a situation applicable to solar coronal con- 
ditions. The plasma motion outside the dissipative layer is de- 
scribed by the set of linear, ideal MHD equations. The wave mo- 
tion inside the dissipative layer is governed by Eq. (ISZb . This 
equation is linear, despite taking into consideration (quadratic 
and cubic) nonlinearity. The Hall terms of the induction equa- 
tion in the perpendicular direction relative to the ambient mag- 
netic field cancel each other out. 

The nonlinear corrections were calculated to explain why 
Eq. (ISZb . describing the nonlinear behaviour of wave dynamics, 
is always linear. We found that, in the second order of approxi- 
mation, magnetoacoustic modes are excited by the perturbations 
of the linear order of approximation. These secondary waves act 
to counteract the small pressure and density variations created 
by the first order terms. In addition, these waves are not resonant 
in the Alfven dissipative layer. In the third order approximation 
the perturbations become Alfvenic, however, these perturbations 



are much smaller than those in the linear order of approxima- 
tion. Equations ( 11081 ) and ( |109b describe the expansion of large 
and small variables, respectively, and demonstrate that all higher 
order approximations of both large and small variables at the 
Alfven resonance are smaller than the linear order approxima- 
tion, provided condition ( l66b is satisfied. This condition ensures 
that the oscillation amplitude remains small inside the dissipa- 
tive layer. From a naive point of view the linear theory is appli- 
cable as soon as the oscillation amplitude is small. The example 
of slow resonant waves clearly shows that this is not the case. 
The nonlinear effects become important in the slow dissipative 
layer as soon as e ~ R^^^^, i.e. as soon as the oscillation ampli- 
tude in the dissipative layer, which is of the order of eR^^^, is of 
the order of '^^ <k 1. We also found that any dispersive effect 
due to the consideration of ions' inertial length (Hall effect) is 
absent from the governing equation. 

This calculation of nonlinear corrections to resonant Alfven 
waves in dissipative layers allows us to apply the akeady well- 
known linear theory for studying resonant Alfven waves in the 
solar corona with great accuracy, where the governing equa- 
tion, jump conditions and the absorption of wave energy are 
alrea dy derived (se e e.g., ISakurai et al.l 119911; iGoossens et al.l 
119951: lErdelvilll998h . 

It is interesting to note that this work can be transferred 
to isotropic plasma rather easily. Shear viscosity, supplied by 
Braginskii's viscosity tensor (see Appendix A), acts exactly as 
isotropic viscosity. Therefore, replacing rj by po^v in Eq. ( l52b 
provides the required governing equation for resonant Alfven 
waves in isotropic plasmas. Moreover, the work on the non- 
linear corrections presented in this paper is also unaltered by 
anisotropy. This implies that we can consider resonant Alfven 
waves in dissipative layers throughout the solar atmosphere and 
still use linear theory if condition ( 1661 ) is satisfied. 
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Technology Facilities Council) for the financial support provided. I.B. was finan- 
cially supported by NFS Hungary (OTKA, K67746) and The National University 
Research Council Romania (CNCSIS-PN-II/53 1/2007) 

Appendix A: Braginskii's viscosity tensor and 
derivation of largest terms 

In this Appendix, we shall derive the largest terms of 
Braginskii's viscosity tensor inside the Alfven dissipative layer 
to be used to study the nonlinearity effects. Braginskii's viscos- 
ity tensor comprises of five terms. Its divergence can be written 
as (iBraginskill965l) 

V-S = %V-So+^iV-Si -h772V-S2-^3V-S3-^4V-S4, (Al) 

Note that the terms proportional to rj^, rj^ and 1)2 in Eq. dAlb 
describe viscous dissipation, while terms proportional to 773 and 
^4 are non-dissipative and describe the wave dispersion related 
to the finite ion gyroradius. They will be ignored in what follows. 
For simplicity in the body of the paper, we have taken 77 = 771 . 
The quantities So, Si and S2 are given by 

So = ^b' (g) b' - ^/j [3b' ■ V(b' ■ v) - V ■ v] , (A2) 

Si = V (gi V + (V » v)^ - b' ® W - W ® b' 

+ (b' ® b' - /) V ■ V + (b' ® b' -I- /) b' ■ V (b' ■ v) , (A3) 

S2 = b' ® W -H W ® b' - 4 (b' ® b') b' ■ V (b' ■ v) , (A4) 
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W = V(b'-v) + (b'-V)v. 



(A5) 



Here v = (u, v, w) is the velocity, b' - Rq/Bq, I is the unit tensor 
and ® indicates the dyadic product of two vectors. The super- 
script 'T' denotes a transposed tensor. 

The first viscosity coefficient, 7]q, (compressional vis- 
cosity) has the follo wing approximate expression (see e.g., 
iRuderman et"ani2000h 



pokBToTi 
% = — 



(A6) 



where po and To are the equilibrium density and pressure, nip is 
the proton mass, kg the Boltzmann constant and t, the ion colli- 
sion time. The other viscosity coefficients depend on the quantity 
w,T,, where w, it the ion gyrofrequency. When w/t, » 1 these 
coefficients are given by the approximate expressions 



1i = 



4(w,T,)- 



'72 = 47/1. 



(A7) 



The viscosity described by the sum of the second and third 
terms in Eq. ( lAll ) is the shear viscosity. For typical coronal con- 
ditions w,T, is of the order of 10^ - 10^, so according to Eq. (|A7| i 
the term proportional to tjq in Eq. ( lAll i is much larger than the 
second and third terms. However, it has been long understood 
that the compression al viscosity does not remove the Alfven 
singularity (see e.g., lErdelvi & Goossenslll995l; iMocanu et alj 
[2008) while shear viscosity does. 

First, we shall calculate the components of the compres- 
sional viscosity. We will use the notation of parallel and perpen- 
dicular components as defined in the paper It is straightforward 
to obtain that 

%(V-S„), = 0, (A8) 



j7o(V-So)x =0, 
?7o (V ■ So)|| = /?o cos a 1 cos a ■ 



(A9) 



c?2 



v± 



■ smff 



dxdz dz?- 

(AlO) 

The shear viscosity, as stated above, is the sum of the second 
and third terms of Eq. ( lAll ). To evaluate these terms we use the 
approximate expression for rj^ and 1)2 given by Eq. (IA7b . As a 
result we obtain 



T), [(V-Si), + 4(V-S2)J=^i 



+(1 -t- Scos^a)— r- +4 



dz^ dxdz 



(All) 



^1 [(V -81)^+4(7 -82) J = 
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(A12) 
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V / /J72 



dxdz 



(A13) 



Equations ( |A8I i-( |AT3T i are complicated, but we can simplify 
them further by taking the largest term only in each equation. 
For the viscosity in the parallel direction it would, at first, seem 
obvious that the largest term will be proportional to \) rather 
than 771 . However, some of the variables proportional to 771 have 
derivatives with respect to x which produce enormous gradients 
in the dissipative layer when there is a transversal inhomogene- 
ity, hence some of the terms proportional to 77; are of the same 
order as or larger than the terms proportional to ^g- It is also 
important to note that in the first order approximation the sec- 
ond and third terms on the right-hand side of Eq. (lAlOb cancel 
(see Eq. |49] l. For the normal and perpendicular components of 
viscosity, the treatment is slightly simpler The compressional 
viscosity is zero, and as derivatives with respect to z are much 
small that derivatives with respect to x, we can select the largest 
term proportional to by observation. Therefore, the viscosity 
tensor can be approximated by 



d\\ 



(A14) 



(A15) 



(A16) 



Equations ( IA14b -( IA16l l give an appropriate approximation to 
Braginskii's viscosity tensor when studying nonlinear resonant 
Alfven waves in dissipative layers. It is interesting to note that 
the terms in Eqs. ( IA14b -( fAT6] l are identical to the largest terms 
when considering isotropic viscosity. Obviously, compressional 
viscosity cannot remove the Alfven singularity since Eq. ( |A9b is 
identically zero. 



Appendix B: The derivation of the Hall term in the 
induction equation for Alfven resonant waves 

In this Appendix we will derive the components of the Hall term 
in the induction equation and show that neglecting the Hall ef- 
fect at the Alfven resonance is a good approximation for typical 
conditions throughout the solar atmosphere. The main reasons 
qualitatively are as follows. When we are in the lower solar at- 
mosphere (e.g. solar photosphere) the Hall conduction is much 
smaller than the direct conduction since the product of the elec- 
tron gyrofrequency, w^, and collision time, t^, is less than unity 
(see e.g.. Priest 1984). For the upper atmosphere (e.g. chromo- 
sphere, corona), where the product cij^Te is greater than unity, the 
Hall conduction has to be considered. However, when the Hall 
terms are derived, the largest terms in the perpendicular direc- 
tion relative to the ambient magnetic field cancel leaving only 
higher order approximation terms which are far smaller than the 
direct conduction. As the dominant dynamics of resonant Alfven 
waves in dissipative layer resides in the components of veloc- 
ity and magnetic field perturbation in the perpendicular direc- 
tion relative to the background magnetic field we can neglect the 
Hall conduction completely from the analysis without affecting 
the governing equation. 

In order to estimate the relative importance of the Hall term 
and resistive term in the dissipative lay er we follow the sophisti- 
cated analysis which was presented bv iRuderman et all ( Il997bl) 
and Clack & Ballai (2008). We do not write down all the steps 
of the analysis, but rather give the salient points specific to the 
Hall effect at the Alfven resonance. 
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Equations ( |29] l and ( l42l i provide the following estimations in 
the dissipative layer: 



,1/2 



) h 



(Bl) 



where e still denotes the dimensionless amplitude of oscillations 
far away from the dissipative layer The thickness of the dissipa- 
tive layer divided by the characteristic scale of inhomogeneity is 
6clknh - This gives rise to the estimations 



ax oz 



th^-Q^-^ih), (B2) 



where h denotes any of the quantities m, bx, b\\, bj_ or v^. Since 
the first term in the expansion of vy is independent of x, it follows 
that 



hnh^ = ^(V||), 



OZ 



(B3) 

We now need to calculate the components of the vectors of 
the resistive term and the Hall term from the induction equa- 
tion normal to the magnetic surfaces (the x-direction) and in the 
magnetic surfaces parallel and perpendicular to the equilibrium 
magnetic field lines. We use Eqs. (IB2b and ( IB3b in order to esti- 
mate all the terms and we only retain the largest. As a result we 
have 



d^b^ 



5x2 



- , -52/7,, 



B() cos^ a d^bj_ 
Bo I 1 dBodb_, 



fioerie \ Bq dx dx 
Bq cos a d^b± 
j-toerie dzdx 



+ cos a- 



dzdx 



+ . . 



+ . . 



(B4) 
(B5) 
(B6) 
(B7) 
(B8) 
(B9) 



where the dots indicate terms much smaller than those shown 
expHcitly. With the aid of Eqs. (IBTT i, ( lB2l i and ( |B3T l we obtain 
the ratios 



AV^-B, 
AV^B_L 
AV'-B« 



(BIO) 
(BID 
(B12) 



For the Hall conduction to be significant in the direction of the 
dominant dynamics of resonant Alfven waves (i.e. in the perpen- 
dicular direction) we must have eoJeTe > 1 . This is plausible for 
the solar upper atmosphere. If this condition holds, then we must 
consider the Hall term in the induction equation. However, if Eq. 
( |B8t is expanded using Eq. ( |29] l from Section III we obtain 



3/2 1 I 1 / '^^o\ db\ 
e A 



Bo\dx I 



(1) 



2a(1) 



d'b 



+ cos a- 



ded^ 



+ ^(eO. 



(B13) 



It should be noted that in deriving Eq. ( IB 13b we have used the 
assumption that ea>eTe - ^(1). The terms inside the braces are 
of the same order as the direct conduction. Hence, they would be 
expected to appear in the governing equation. When substituting 
for b^;P and b^^^^ using Eqs. (05]) and ( l46b . respectively, it is found 
that the terms inside the brackets cancel 



cos a I dB()\ 5m*'' cos a I dB(A 5m*'' 



V \dx j V \dx I 



= 0. 



(B14) 



Equation (IB 14b shows that the Hall term in the perpen- 
dicular component of induction is always smaller than the di- 
rect conduction in the solar atmosphere. The normal and paral- 
lel components of the Hall conduction are, in fact, larger than 
the perpendicular component. Nevertheless they play no role in 
derivation of the governing equation of resonant Alfven waves 
in dissipative layer. The parallel component is the largest of the 
three components and this is to be expected as the Hall eff'ect is 
strongest at right angles to the dominant wave motion. This is in 
complete agreement w ith the study on resonant slow waves by 
IClack & Ballail (|2008|) which found the largest Hall effect was 
in the perpendicular component of Hall conduction, which is at 
right angles to the dominant dynamics of resonant slow waves. 

In summary, it is a good approximation to neglect the Hall 
term in the induction equation when studying resonant Alfven 
waves in dissipative layer. This approximation holds throughout 
the entire solar atmosphere. 
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